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Математическое моделирование процессов, 
разрьівньх на линиях триангуляции 


У статті побудована математична модель розривного процесу, що описується функцією двох змінних, 
у вигляді розривного інтерполяційного лінійного сплайна. Причому функція має можливі розриви на 
лініях триангуляції довільними трикутниками. Побудовані розривні конструкції включають в себе, як 
частинний випадок, класичні неперервні сплайни. 

Ключові слова: розривний процес, розривна інтерполяція, 

триангуляція довільними трикутниками. 
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В статье построена математическая модель разрьтвного процесса, которьтй описьтваєтся функцией двух 
переменньхх, в виде разрьівного интерполяционного линейного сплайна. Причем функция имеет возможньте 
разрьвьт на линиях триангуляции произвольньми треугольниками. Построеннье разрьївнье конструкции 
включают в себя, как частньй случай, классические непрерьівнь6е сплайнь. 

Ключевьге слова: разрьвньгй процессе, разрьівная интерполяция, 

триангуляция произвольньми треугольниками. 


Вступ 


Задача наближення розривних функцій є однією з найскладніших задач обчислю- 
вальної математики. Спеціалістам з обчислювальної математики добре відомі оператори 
наближення неперервних та диференційовних функцій за допомогою поліномів та 
сплайнів (1), 12). Відомі також праці з наближення неперервних функцій однієї змінної 
кусково-сталими функціями |31, 4), в яких неперервні та диференційовані функції 
наближуються сплайнами степеня нуль. В роботі |5| була розглянута апроксимація 
розривних розв'язків (функцій однієї змінної) диференціальних рівнянь за допомогою 
розривного методу Гальоркіна. А в роботі |6) розглядається розривний метод Гальоркіна 
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для еліптичної крайової задачі з використанням двовимірних неузгоджених сіток. 
Цей метод дозволяє враховувати неконформність елементів. Причому метод забезпечує 
неперервність розв'язку, хоча від базисних функцій узгодженості не вимагає. 

Таким чином, у вказаних роботах досліджувалося наближення неперервних функ- 
цій за допомогою неперервних та розривних сплайнів або розривних функцій за допо- 
могою неперервних. Але загальної теорії таких наближень не існує. В даній роботі 
ми пропонуємо таку загальну теорію побудови розривних сплайнів, множина яких як 
частинний випадок включає множину неперервних сплайнів, що можуть мати розриви 
першого роду в заданих точках або на заданій множині ліній - границь елементів. 

Задачі наближення розривних функцій виникають частіше, ніж задачі наближення 
неперервних функцій. Наприклад, в методах комп'ютерної томографії на даний час 
недостатньо вивчене питання про використання інформації про внутрішню структуру 
тіла людини (різні органи мають свою форму та щільність тканин). 

Тобто актуальною є розробка та дослідження теорії наближення розривних функ- 
цій за допомогою розривних функцій. 

Авторами вже досліджувалися питання наближення розривних функцій роз- 
ривними сплайнами з використанням елементів, що мають хоча б один прямий кут. 
В роботі |7| був запропонований метод наближення розривних функцій двох змінних 
розривними інтерполяційними білінійними сплайнами, а в роботі |8| - інтерлінацій- 
ними розривними сплайнами на ректангульованій області визначення. Були також 
побудовані розривні інтерлінаційні сплайни для наближення функцій двох змінних, 
область визначення яких розбивається на прямокутні трикутники |З). 

Дана стаття присвячена побудові та дослідженню інтерполяційних розривних 
сплайнів для наближення розривних функцій з областю визначення, що триангулюється 
довільними трикутниками. 

Метою даної роботи є модифікація запропонованого методу сегментації символів 
у напрямку збільшення його якісних показників в умовах перспективних спотворень 
зображень номерних знаків та їх часткового затінення. 


Постановка задачі 


Нехай в області р - 017 розривний процес описується розривною функцією 
ЛО, у). Нехай Р розбита на п довільних трикутників, та інформація про функцію 
задана у вигляді значень функції /(х, у) у вузлах трикутників області І). Функція 
ДО, у) має розриви першого роду на границях між цими трикутниками (не обов'язково 


між всіма). Метою роботи є побудова та дослідження математичної моделі розривного 
процесу, що описується розривною функцією /(Х, у), у вигляді оператора розривної 
лінійної сплайн-інтерполяції. 


Побудова математичної моделі розривного процесу 


Розглянемо довільний трикутник Т., ізіп з вузлами А (х9, у, і 21,п, 


Кк - 1,3 (рис. 1). Вважаємо, що на кожній із сторін заданого трикутника функція /(Х, у) 
може мати (а може і не мати) розриви першого роду, причому у вершинах трикутника 
функція набуває значень 


СО - | то | Р(х,у), ізіп, б -1,2,3:- 
(х, у) А 
(х,у)єтТ; 
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Рисунок 1 - Зображення довільного трикутника з триангульованої області 
визначення розривної функції /(Х, у) 


Визначення. Будемо називати розривним інтерполяційним лінійним поліноміаль- 
ним сплайном в області Т. с 0), і -І,п наступну функцію 


0) 0) (1 
5; С, зеніт (2) 9 (5 У) | о 9; (о У) 


об) (А?) і о (А?) Р (1) 


де 
ка) (1 
хоХ; яз У; 
яко) (ху)ву- уФ | ; 209) У 
- - 2-х 
і і 
30) (1 
хх; зсуву, 
Фо (ху) у- у о МЕМ аа ем, т й, - 
Х; бе Х; 


(2) (3) (2) 
3) з зу Кох 02 -У) 
Фо РОдре НУ о о 
і і х9) М хО) 
Теорема 1. Функція 5(х, у) - 5, (Х, У), (х, у) є Т, с ДЮ задовольняє інтерполяцій- 
ним властивостям 
т 5. (х, у є С, іа1,п, К -1,2,3. 


(зу зАМ 
Со у)єт; 


Доведення. Перевіримо твердження теореми для К- 1. 


Ши 5 (х,у)- Ши 5,(х, у) - 
(ху) А? зх? 
(х, у)єТ; узуФ 


Зх УЧ УМО ЕН СРО РО 
осо ООН) | су ФОН) | сі ФРС У 


З ЗА є В ри о В 
- 0) (9 ; (2) (2) ; (9 0) 
(АР) о (А?) (47) 
3, а) (а 223 Й ПРОЛНІ Є Б ОО 
і) Фо (39) (2) о е (3) они баб мо МИ 


о АФ?) о? 40) ЕС У 40)) 
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І І 3 І 
есрос ФО УФ. с ох у) І б 
і і і і хо) Ми жу і?) (А?) 
І І 7 І 
ес Фо хі) зх - у) 1 осі. 
і ПОР 9 а) (А2)) о 


Теорема доведена. 
Наведемо оцінку похибки наближення функції лінійним інтерполяційним сплай- 
ном, яка наведена в роботі Суботіна (10). 


Нехай совОвУреНУ, вро Ре сі, р Сх, У) 2 6 аа, ше 


похідна за напрямком 8. 
Теорема 2. Нехай /(х, у) є М., Мо «(ЛО У Де Г (о у) - неперервні в Д та 


ра) В.Г є 


ратором 5(х, у) - 5; (Х, у), (х, у) ЄТ; с Р, тоді для оцінки похибки наближення в кож- 


, Ми з(циз) є ), Ху « (у, у,) є ), М5) наближується опе- 


ному трикутному елементі розбиття справедлива они 
(РО, У) З 505, У) є «Мій, 


де п - найбільша із сторін трикутника. 
Теорема 3. Якщо С? - Г(А9), Кк «13,і-1 п, то в кожному трикутнику Т,, 


і - І, п, оператор (1) точно відновлює всі лінійні функції. 

Доведення випливає з того, що через три точки можна провести тільки одну 
площину. 

Зауваження 1. Якщо значення функції у вузлах одного трикутника збігаються зі 
значеннями функції в цих же вузлах зі сторони сусідніх трикутників, то оператор (1) 
є класичним неперервним лінійним інтерполяційним сплайном. 

Зауваження 2. Якщо значення функції у вузлах трикутної сітки невідомі, то для 
знаходження невідомих коефіцієнтів С і ), Ка12,3 і-іп вданій роботі пропону- 
ється використовувати метод найменших квадратів, згідно з яким всі невідомі зна- 
ходяться з умови 

2 . 
ЛІСУ ь3 ПІ /0ь у) -з, (ху, С) | хау -» піт. (2) 
тер т; 

І тоді отримаємо апроксимаційний розривний лінійний сплайн. 

Оператор (1) з коефіцієнтами, визначеними формулою (2), являє собою матема- 
тичну модель процесу, що має розриви на лініях триангуляції області визначення 
функції двох змінних, яка і описує згаданий процес. 


Приклад. Нехай задані вузли триангуляції одиничного квадрата Р - |0; Ще», 


0 0 0 1 0.3 | 
0 1 03 ол 1 1 
Ро , У1з НО зраднОО уро ода , УЗа ї 
0.6 04 1 1 1 0 
0 0 0.6 04 1 0 
0 0 03 ол 0.6 04 
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Тобто одиничний квадрат розбито на 6 довільних трикутників (рис. 2а) 
7 
ТІ аоерюєхкоз Пааукічя| 


Т. «|соера-уєхкозі зо-0тх олкуєі) З 
3 
Т. ерюзаановазяє у«ольбаодо 
3 6; 
м(х,уує р:0,6«х«ї; КН РОН УВА ОК е 


Ту «|серюбехамсяє у«оачбаноє| : 


п (аудерідуєачі-уочуєолі) 
2 3 2 
п одерюєхсозакеуєміо аорерозехаолдакуєі-х), 


Нехай в області Ї) задана розривна функція /(х, у), яка має розриви на лініях 
заданої триангуляції, але не на всіх (рис. 26). 


І, окха0б дякуєі-х 
їі 
Х, 0,7«у«Ї чав юю 


ДО» у)є ; . 
3-х, ЧОРНО ОРОНИГО ВРО 


-2х'42, осу«0А бука 


Функція /(х, у) має розриви першого роду у вузлах заданої трикутної сітки та 
в них має такі значення: 


пт о (ху) -Ї; йт  Г(х,у)-Ї пт /(х,у) 2; 
(5 У)- (0,0) (5. У)- (0,0) (5 У)-2(0,0) 
Со У)єТ, (5 У)ЄТо Со у)єТ; 
йт  /(х2,у)-Ї йт  /(х2,у)-0); йт ЛО, у) єї 
Соу)-м(О) (ху) -м(О,1) (Со у)-М0.3,0.7) 
Со У)єТ, (5 У)єТ, Со У)єТ, 
т Т(х, у) - 0.3; йт Тео ууАТ йт ГО у)-б 
(Со у)-(0.3,0.7) (Со у)-0.3,0.1) (Со у)-М0.3,0.1) 
Со у)єТ» Со у)єТз (о У)ЄТо 
йт Рогу) В йт ЛО, у) 22.4 йт Г(х, у) -2.4; 
(ху) М0.6,0.4) (Со у)-м(0.6,0.4) (ху) М(0.4.0.6) 
Со у)єТо С у)єТз Со у)єТ, 
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Пт ЛО, у) - 1,65; пт /(х,у)-Г пт /(х,у)- 2; 
(х, у) (0.4,0.6) (ху) (ху) 1) 
С у)єТ; С у)єТ, Со у)єТз 
Жов уРеД пт 7/(х,у)- 23; йто /(х,у)-0. 
(5 У)- (11) С у) Ц,0) (5 у) 0) 
Со у)єТ; Со уує Ту Сх, у)єТ; 


Ху 


Рисунок 2 - а) триангуляція області визначення функції /(х, у); б) графічний 
вигляд наближуваної функції /(х, у) 
В кожному трикутному елементі побудуємо інтерполяційний сплайн 5(х, у) У 


вигляді формули (1). Після знаходження за формулою (2) невідомих коефіцієнтів, отри- 
маємо сплайн вигляду (рис. 3). 


1, окха06 Закуєічя 
7 
З 0,7«ук«Ї 09840600 
5(5, у) - 3 
3-х, оо п 


-2.1х-0,.5у-24, 0«у« 0.4 5 «х«І-у 


55 


Рисунок 3 - Побудований наближуючий розривний сплайн 
Побудований розривний сплайн точно наближує ту частину функції, де вона є 
постійною або лінійною, що і підтверджує викладену вище теорію. 
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Висновки 


У роботі запропонована математична модель розривного процесу у вигляді роз- 
ривного інтерполяційного лінійного сплайну. Вважається, що розривний процес описує- 
ться функцією двох змінних з розривами першого роду на лініях триангуляції довільними 
трикутниками. Метод, побудований авторами, припускає, що розриви наближуваної 
функції відомі, і тому вони збігаються з розривами наближуючого сплайна. Наступним 
кроком планується розробити методи наближення розривних функцій розривними 
сплайнами, коли розриви наближуваної функції ще треба знайти. А також планується 
застосувати розроблену теорію наближення розривних функцій розривними сплайнами 
до розв'язання двовимірної задачі комп'ютерної томографії. 
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КЕЗ5ЗСМЕ 


О.М. Іуїууп, У.І. Рег5Ппіпа 
Матетайса! Моаєгіпе ої Ргосе55е5, 


Різсопипиодиз оп Тгіапеціатоп І іпез 

І агісіе Ше ргобіет ої соп5ігисйоп ої пзаїретайсаї піодеї ої Фізсопіпиоиц5 
ргосе85 мрісП 15 дезсгібед Бу Фіізсопійписиц5 Гипсоп ої (мо уагіабіе5 15 соп5ідегеай. І 15 
зиррозей, Шас Фе Гипсйоп гапее ої дебпійоп БгеаК5 іпіо апу ііапеЇе5 (пої песеб5агіЇу 
гесіапешаг) уПпісп аге пої риє еасп оїег, апа Шеїг рагіез аге пої сго85ед. АЇЗ8О їй 15 
зирроз8ей, Факс Типсйоп Паз8 роз58ібіе гиріиге8 ої Ше Пг5і 5огі оп Пппез Бебуееп (гіапеціаг 
еіетепіз. 

Тре таШетайса! поаеі 15 ипдег соп5ігисйоп їп ре Їїогт ОЇ йзсопипиоц5 Ппеаг а 
зріпе-іпіероїацоп. І ууогк Ше е5ійтайшоп ої Ше соп5ігисіва фізсопійпиоця дезієп їп еасп 
ішіапошаг еР|етепі 15 ге5ийей. Пп у/огк сопдаійопя ас мрбіср іпіегроіайопа! дфізсопипиои5 
зріпе 15 сопійпиоц5я аге гезиіей, Па! 15 Ше єепега! Шеогу ої арргоасП ої Ше 5еї Гипсйоп5 
міс 5ресіа! са5е 15 а Ше Кпомп Шеогу ої арргоасі Бу сіазбіса! сопіїпиоця 5ріпез 15 
деуеіорей. ТРе ехатріе ої арргоасі ої й4їзсопіпиоця Гипсйоп ої (м/о уапіабіе5 Бу соп8ігисівд 
іпегроіацопаї Фізсопіпиоцз єріте 15 гезиїей, апа Ше гесетуеа гезиів ої арргоасі аге апа|у7гед. 

Айфог58 5пом, Шаї Ффзсопійпиоця їп 5оте роїпі5 ог оп 5оте Ппе58 ої Гипсйоп оп 
сумо уагіабіез її 15 Бейег (о арргоасп Фізсопіпиойя зріпез. ТПиє, 16 15 ро85іБІе (0 гесеїуе 
едиаПу арргесіайоп ої ап еггог ої арргоасі іп еасП еЇетепі ої зрійцтя, мурісі аге іпрегепі 
їп 15 сопіпиоц5-діНегеппатеа з5ріпез. 


Стаття надійшла до редакції 01.06.2012. 
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